Theory of spin dynamics of magnetic adatoms traced by time-resolved 

scanning tunneling spectroscopy 

M. Schiller, Y. Pavlyukh, and J. Berakdar* 

'^' I Institut fur Physik, Martin-Luther Universitdt Halle Wittenberg, 

O ' Halle, Heinrich-Damerow-Strafie 4, 06120, Germany'^ 

(N 
^ . (Dated: March 1,2013) 

Abstract 

pi . The inelastic scanning tunneling microscopy (STM) has been shown recently (Loth et al. Science 329, 



7j 



1628 (2010)) to be extendable as to access the nanosecond, spin-resolved dynamics of magnetic adatoms and 
molecules. Here we analyze theoretically this novel tool by considering the time-resolved spin dynamics of 
a single adsorbed Fe atom excited by a tunneling current pulse from a spin-polarized STM tip. The adatom 
j^ . spin-configuration can be controlled and probed by applying voltage pulses between the substrate and the 

fi \ spin-polarized STM tip. We demonstrate how, in a pump-probe manner, the relaxation dynamics of the 



sample spin is manifested in the spin-dependent tunneling current. Our model calculations are based on the 



C/3 

^ ■ 

p . scattering theory in a wave-packet formulation. The scheme is nonpertubative and hence, is valid for all 



voltages. The numerical results for the tunneling probability and the conductance are contrasted with the 
^ ' prediction of simple analytical models and compared with experiments. 
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I. INTRODUCTION 

Spin systems are extensively utilized as essential elements for (quantum) information storage or 
computing devices [1,2]. Thereby, environmental effects play an important role as they generally 
lead to decoherence and relaxation. The time scales for these processes depend on the underlying 
coupling mechanisms and exhibit in general a marked temperature and dimensionally dependence. 
A detailed insight and a possible control of this relaxation is a key factor for the operation of these 
devices. Desirably the relaxation time should be larger than the operation (switching) time. In 
this sense, it is of importance to identify spin systems with suitable relaxation properties and 
amenable to full fast control of the spins. Prototypical examples are realized in semiconductor 
nanostructures [3, 4] or nanomagnets [5, 6]. For a single magnetic atom or molecules adsorbed on 
a substrate, the surface spin excitations are usually weakly coupled to bulk magnons resulting in 
relaxation times on the order of nanoseconds, which has been demonstrated to be sufficient for a 
coherent spin manipulation [7]. 

The experimental analysis of surface-deposited structures can be performed conveniently by 
means of the scanning tunneling microscopy, which has an advantage of an atomic spatial res- 



olution. Using a spin-polarized tip, the spin of a single magnetic adatom can be probed and 
manipulated [8], since the tunneling probability and the current depend on the relative alignment 
of the surface and the tip magnetic moments. Very recently the same group demonstrated the po- 
tential of the technique for atomic-scale information storage and retrieval [9]. The time resolved 
STM experiments [10] renders possible the observation of the quantum dynamics. The transient 
precessional dynamics of the excited spin states is still too fast to be measured. The relaxation 
process can, however, be mapped onto the current-dependence in a pump-probe experiment. Thus 
one has a possibility to directly measure the relaxation time of such systems which opens the way 
for testing different configurations with maximal coherence time. 

To our knowledge, the first experiment to measure the relaxation times via STM has been 
carried out by S. Loth et. al. (Ref. 8) for a magnetic system with a particularly long spin relaxation 
time of above 200 ns. The basic experimental setup (figure 1) consists of a Fe-Cu dimer placed on 
the Cu (100) surface covered with a CU2N overlay er and then probed by a spin-polarized tip. A 
magnetic field of i? = 7 T is applied, almost aligning the spin of the Fe atom parallel to the spin 
axis of the tip. Due to the magnetic coupling between the spin of the tunneling electrons and the 
spin of the Fe atom, the adsorbate is driven into an excited state with a different projection of the 
spin moment on the spin axis of the tip. This influences the tunneling current and allows to trace 
the spin dynamics. 

We will consider theoretically the experimental situation (section II) using scattering theory 
(section EI) with the aim to calculate the spin-dependent transport properties. The method is com- 
plementary to other formalisms. A widely spread approach is to describe the tunneling by the 
master equation with transition rates calculated in a standard "golden rule" manner [11, 12]. This 
method is, in essence, a perturbative formalism and hence requires a separate justification for each 
of the considered experimental configurations. The nonequilibrium Green's function approach has 
a potential to systematically take into account many-body effects, however, the transition rates are 
often introduced as adjustable parameters [13]. In the present method we focus on the descrip- 
tion of the tunneling starting from a model Hamiltonian including the magnetic anisotropy, the 
exchange coupling and the coupling to the external magnetic field. Our theory proves that it is 
indeed possible to directly measure the spin-relaxation times provided they are longer than the 
precession time. Furthermore, our main findings for the relative number of transmitted electrons 
with respect to the pump-probe delay (presented in section V) allow for a rigorous determination 
of the magnetic coupling parameters. 
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FIG. 1. Color online (a) A spin-polarized STM tip is placed above a single Fe adatom adsorbed on the copper 
(100) surface with an CU2N overlayer to excite and probe the spin dynamics of the adatom. The tip consists 
of one Mn atom attached at the end of a tungsten tip. The Mn atom acts as a spin filter for the electrons 
propagating from the tip to the sample. A strong magnetic field of B = 7 T adjusts the spin projection of 
the Fe atom nearly parallel to the spin polarization direction of the tip. (b) A theoretical model of (a) by a 
tunneling barrier and a spin-dependent interaction potential around the Fe atom on the surface. The solid 
line corresponds to the surface-deposited ground state spin. The coupling between surface-deposited spin 
and tip-spin states is antiferromagnetic. Hence, the coupling energy decreases when exciting the Fe spin 
states. The effective tunneling barrier for the first excited spin state is indicated by the dashed line. 

n. MODEL HAMILTONIAN 



Our model Hamiltonian matching the experimental situation in figure 1 includes 

2 
2m 



(1) 



which describes the tunneling electrons (with a mass m and momentum operator p) emanating 
from the tip and are subject to the magnetic field B. They are coupled to the surface spin subsystem 
via a spin-dependent term Vb- The spin of the tip electrons is described by the vector operators 
f = {tx, Ty, Tz) containing the Pauli matrices in their standard representation. In a similar way, 
S = {Sx, Sy, Sz) are the operators for the surface spin. The tunneling barrier Vb{x) is measured 
relative to the Fermi levels of the tip or the sample (for simplicity we assume both to be equal). 
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TABLE I. The spectrum of the adsorbate spin Hamiltonian (2), referred to the ground state, and the projec- 
tion of the spin along the magnetic field. The third column contains the quantum number of |m) states with 
the dominant contributions to \xu)', the values are given in the last column. 
1/ S^-Ef(meV) {Xu\Sz\Xu) m \{m\xu)\'^ 

0.855 
0.869 
0.544 
0.544 
0.944 

The last term in equation (1) describes the spin of the tip electrons in the presence of the magnetic 
field B. 

The second subsystem is the surface spin in an anisotropic environment (due to the broken 
translational symmetry and other atoms placed nearby). The appropriate Hamiltonian for this part 
reads [5, 14] 

Hs = g^iBBS, + 5S', + e{Sl - S'^), (2) 

which indicates that the easy axis is parallel to the (strong) magnetic field. We have chosen typical 
values for a rather strong magnetic anisotropy. This does not directly resemble the experimental 
situation, but serves to illustrate the effects arising from the anisotropic environment more clearly. 
In particular, we take 5 ~ — 5 meV and e ^ 1.5 meV for the anisotropy constants. The value 
of the magnetic field is kept with 7 T. The Lande factor g is experimentally [14] found to be 
approximately 2, which corresponds to the total spin of the free Fe atom (S = 2). The Hamiltonian 
(2) has 5 eigenstates \xu) with eigenenergies E^ (table I lists their properties). The ground state 
in orbital representation is depicted in figure 2. 

The third term on the left side in equation (2) results in a mixing of the different eigenstates 
\m) of the free spin in the magnetic field. Thus, the expectation value {S^) in the ground state 
deviates from the noninteracting value -2. The ground state contains a large amount of free spin 
configurations (see table I). The first excited state corresponds to the spin pointing in the opposite 
direction. Hence, the relaxation \x2) -^ \Xi) is dominated by the transitions \m = +2) —^\m = 
-2). 

For the spin-dependent interaction of the tunneling electrons with the adsorbate we adopt the 




FIG. 2. The ground state xi(^) (p) ^^ orbital representation. The shape is given by \xi{^j 0)|> whereas the 
color encodes the real part of the wave function. The red arrow denotes the expectation value of the spin 

(S). 



picture of extended states (tunneling electrons) coupled to a localized magnetic moment. In brief, 
such a process is known to be determined by the Hunds rules and the local electronic correlation. 
E.g. for s states (with energy e^) coupled to a pinned magnetic impurity (with energy ed) one may 
utilize the Anderson model [15]: if And = Efc,^ ^fc^^fca + Ea (^dnda + Und^ridi + Hds, where a labels 
the spins and Uka (nka) is the corresponding particle density operator. U stands for the Coulomb 
repulsion between electrons with antiparallel spins. The term H^s = Tlik a^kdCk^Cda + h-C ac- 
counts for the s-d hybridization with the mixing matrix element Vkd- The Anderson Hamiltonian 
is a model for a delocalized s-band interacting with a single (i-orbital. As we are interested in 
the quantum dynamics of the local magnetic impurity we should map the Anderson model onto 
the Kondo Hamiltonian [16] (-ffKondo = —Jo^s ■ <^d) with the effective Kondo coupling constant 
J. Tunneling electrons in the vicinity of the Fermi level are of relevance for the tunneling. This 
leads to estimate J = 2\Vkpd\'^ ^ (F+m ^ ^- Typical values [17] are in the order of | V^^rfp ^ 2 eV, 
t/ ?^ 8 eV and e^ ?^ 5 eV, meaning that J sa 1 eV. 

In STM experiments the tip-sample distance (we take 5 A) is usually not negligible to the 
typical extent of (i-orbitals in Fe adatoms (which is about 1 A). In the calculations, we use therefore 
an interaction iJint of the Kondo type with a spatial-dependent Kondo coupling 



ifint = UQW{x)f ■ S. 



(3) 



The Mo > is the coupling strength and the function w{x) is a dimensionless form factor defining 
the width of the interaction region. We choose a strongly localized function with a width about 

lA. 



The total Hamiltonain of the system is thus 

H = Ho + Hs + H,^f (4) 

The total effective potential experienced by the tunneling electrons Vb{x) + {xu, t \Hint\Xu, t) is 
thus different for the ground state u = 1 and the first excited state z/ = 2, as sketched in figure 1 (b) 
(the tip spin is f in both cases). The anti-ferromagnetic coupling of the tip spin and the surface 
spin effectively raises the tunneling barrier for the tip electrons, for they have to overcome the 
additional spin-interaction energy. On the other hand, any excitation of the surface spin lowers the 
effective potential barrier (c. f. solid and dashed lines in figure 1 (b)). Therefore the excitation 
Ixi) — 7- 1x2) is favored as long as the applied voltage is higher than the first excitation energy, i. e. 
if eV > E2 — Ef. For a smaller voltage this excitation is not reached energetically. All higher 
excited states also lower the effective potential barrier, but the effect is maximal for \x2)- This is 
reflected in the enhanced excitation probability for the first excited state. 

If we on the other hand consider the case of | J,) -electrons beeing emitted from the tip, the 
effective potential Vb{x) + {Xv^i \H\at\XvA) undergoes a change of the sign in the last term. 
Therefore, {xi^i |-f^int|Xi)i) becomes a negative contribution and thus enhances the tunneling 
current, whereas (x2!i |-^int|X27i) acts as an additional barrier. The other consequence is the 
smaller excitation probability because the transition \xi) — ?• IX2) does not reduce the tunneling 
barrier. Physically, this picture reflects the anti-ferromagnetic coupling again. 

We remark that the existence of a direct transition |xi) — )■ \x2) is governed by the transition 
matrix element (x/^, t \Hint\Xu, t) i^i the surface spin space, which is connected to the anisotropy 
constant e. For the rather large value of £ ^ 5 meV, the direct excitation |xi) — ?• IX2) is the 
dominant process [18]. 

III. SCATTERING FORMULATION 

In a wave-packet scattering picture of the process depicted in figure 1 one identifies an undis- 
turbed part i/frcc = p^/(2m) + gefiBBrz + Hs and a scattering part Hsc = Vb{x) + Hint- The 
unperturbed (asymptotic) scattering states are cast as (± stand for incoming/outgoing waves) 
1$^^) = I ± A;) (g) \xu) ® \t) where \xu) (|t)) is the spin ground state of the surface (tip). | ± A;) is 
the orbital part of the wave with a wave vector k which is determined by the direction and by the 



amplitude of the applied bias voltage V, i.e. k = y^2meV/h? (and e is the electron charge). By 



propagating |$^^) the state of the eigenstates |^^^) of the full Hamiltonian can be obtained (up 
to a phase) [19-21], i.e. |V'+,) = f/(0, -oo)|<l>+,) and |^-,) = f/(+cx), 0)1$-,), where U{t,to) 
is the propagator in the interaction representation. Tunneling is governed by the probability am- 
plitude to scatter from the incoming eigenstate Itp^a) i^^to the outgoing state \ip^^j.) i^ given by 
the scalar product {'ipu,T\''PXa-) — {^u,TWi+oo,—oo)\^^^^). The last equation defines the on-shell 
S-matrix S = U{+oo, — oo). Evaluating the S-matrix element 

5;: = {^-^\s\^;j (5) 

we find the tunneling probability to scatter from the tip electron with a spin state \t) into \a) 
and from the surface spin state \xu) into \x^) by taking the absolute value of the S-matrix (equa- 
tion (5)). For a given initial spin-spin density matrix po (of the dimension 10): 

we obtain the total tunneling probability T{V) and, thus, the normalized density matrix elements 
by 

p;: = ^($;,,|t/(+oo,0)pof/(0,-oo)|<|.+J. (6) 

The probability to scatter into a specific channel {v, a) is given as 

T,^{V) = ($-,|f/(+oo,0)pof/(0,-oo)|$+,), (7) 

Whereas the total tunneling probability is r(V^) = ^^ , T^^r (^) = Trp. 

In the setup figure 1 we are interested in the surface spin dynamics only. Hence, we trace out 
the spin \a) of the tunneling tip electron. This yields the reduced density matrix of the surface 
spin. Furthermore, we fix the initial spin state of the tip electron as |r) = 1 1) as determined by the 
polarization of the tip and \xu) = Ixi)' i- £■ the surface spin is in the ground state when we switch 
on the voltage. We simplify the notations and write the resulting density matrix as p'^^(V). The 
population of the final surface spin state \xu) as a function of the applied voltage is then given by 
Pui^) = Puui^)^ whereas the phase information can be extracted from the off-diagonal elements. 

The assumption |t) = | f) corresponds to the case of a full polarization of the tip. A partial 
polarization /3 can be included in the formalism via the initial spin density matrix of the tip, i.e. 



For clarity however, we employ (3 = 1 throughout this work and comment on the change of the 
results for /3 < 1. 

To calculate the individual S-matrix elements we utilize the wave packet approach [22, 23]. 
This method exploits the fact that in the Schrodinger picture the full eigenstates can be expressed 
as a superposition of propagated wave packets, i.e. 

where the normalization factor 77^^ is defined by 



The wave vector k for 77+^, i. e. for the incoming waves, is given by A; = ^/2meV/W, whereas for 
the outgoing state, the wave vector may be reduced due to inelastic tunneling. In this case, we set 
k = \ 2m{eV + Ef + Er^ — E^ — Ea^)/h? if the term under the square root is positive. Here 
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-E^'P {E^) is the energy of the tip (surface) with a spin state 7. 

The wave packets |0j^^) are superpositions of full eigenstates. However, by placing the wave 
packets far away from the interaction region ?7(t, to) becomes the unit operator in which case the 
packets by composed of the eigenstates of -fffree, i- e. we can write 10^^,-) as the product state 
\g±)\Xv)\T~) with a Gaussian function g±. The subscripting ± now corresponds to centering the 
Gaussian to the left and to the right of the interaction region. 

Using the spectral representation (equation 9) we find the S-matrix elements as 

The formula (11) is then evaluated by (i) propagating the initial product state, (ii) calculating 
the overlap with the final state on the right side of the barrier, and (iii) by performing a Fourier 
transformation. 

IV. SPIN STATE POPULATION 

We now proceed with computing the population P^ of the surface spin eigenstates \xt) as a 
function of the applied voltage V for different coupling constants uq. We assume the surface spin to 
be initially in the ground state and the spin of the tip's electron is being aligned with the magnetic 
field. Figure 3 depicts the dependence of the population of all surface spin states. The energy 



of the tunneling electrons (eV) determines which inelastic scattering processes are possible. For 
voltages below the excitation threshold the corresponding spin state can not be excited due to the 
energy conservation. For the other transitions, one has to take the spin of the tip's electron into 
account, as well. By evaluating the spin-coupling f ■ S* in the basis of product states one finds 
that the transitions |xi) -> Ixs) ^^^ \Xi) ~^ IXi) ^e only allowed if accompanied by spin flip 
processes. Other transitions require the conservation of the spin direction. We can, thus, define 
specific voltages K, where new tunneling channels open, i. e. eV^ = E^ + E^^^ — E1^^ — Ef. For 
1/ ^ 3, 4 we obtain eK = E^ - Ef whereas for i/ = 3, 4 eK = ^f + E^^ - E"p - Ef. The 
voltages Vy are shown in figure 3 by vertical dashed lines. 

The basic trend of the voltage-dependence of the population is to reduce the coupling energy 
by exciting the surface spin. As long as the energy of the tunneling electrons is high enough, the 
population of the ground state decreases. For voltages in the interval between the first two dashed 
lines, the first excited state is the only accessible tunneling channel, so that the population is solely 
transferred to \x2)- The efficiency of this transition is determined by the coupling constant uq. 
If the voltage increases further, the next tunneling channels open and the spin interaction energy 
can again be lowered by exciting the surface spin. Therefore, P2 decreases and P^ rises by the 
same amount. The situation is similar between the last two dashed lines, where P^ is decreased 
in order to increase P/. The highest state \x^) corresponds to a very small projection parallel to 
the magnetic field as the major contribution is \m = 0) (see table I). The transition probability is, 
thus, much smaller than a transition into the other states. For the surface spin it is therefore more 
convenient to maintain a higher population of IX4) instead of decreasing Pf when the voltage 
exceeds the last threshold. 

For larger coupling constants uq the probability to scatter to the excited states of the surface 
spin increases. However, this holds true as long as uq < Mq, where Uq is about 7 eV. For larger uq, 
the transition efficiency decreases again. In order to reveal the dependence on uq, we calculated 
the population of the ground state as a function of V and Uq (figure 4). Since the effects are 
especially pronounced for the transition \xi) — ?■ IX2), we will only consider the corresponding 
voltage interval. As it turns out, the excitation of the surface spin is only probable for certain 
values of the coupling constant. For large uq, the excitation probability even approaches zero. 

In order to explain this behavior, we have to apply a further simplification that renders possible 

10 



an analytical treatment (see appendix). Analyzing the simplified model we infer that 

2mxoC22 

with C22 = {X2 1 'S'2 1 X2 ) and q = ^/2mVBjh? for the optimal value of the coupling strength. Phys- 
ically, the occurrence of this optimal value can simply be interpreted as a resonance phenomenon 
of scattering on the narrow quantum well of the spin-spin-interaction potential. The value of Uq is 
represented by the horizontal white dashed line in figure 4. 

Generally, for values for the tip polarization (3 < I, the excitation probability decreases. We 
have already noted that the transitions induced by the tunneling electrons with the opposite spin 
direction (|r) = ||)) are weaker in comparison to 1 1) -electrons. For (3 decreasing from (3 = 1, the 
strong (It)) and the weak (| i)) contribution are mixed together [24]. Despite from beeing reduced 
to a smaller intervall of the probability, the shape of the curves in figure 3 remains the same for 
realistic values of /3 ^ 30 %. 

The resonance phenomenon arising from the quantum well of the effective potential becomes 
also less visible for a smaller tip polarization, as the quantum well is inverted to a barrier for 
electron spins of the opposite direction. The value of Uq is however hardly influenced. 

V. RELAXATION DYNAMICS 

STM experiments are not only suitable for determining the transport properties, but also offers 
an insight in the relaxation dynamics of the excited surface spin by using a pump-probe scheme. 
The time resolution is not sufficient to resolve the precession or the transitions between neighbor- 
ing |m) states, but it is still high enough to measure the relaxation from \m = +2) to \m = —2). 
This corresponds in a good approximation to the transition IX2) — ^ |xi)- 

The idea of determining the relaxation dynamics is to apply a pump voltage pulse in order to 
excite the surface spin, i. e. to increase the population P2 and to detect the time evolution by a 
second probe pulse. The voltage of the latter has to be chosen according to the excitation threshold. 
We set the maximal values Vpump = 3 mV and V^^.^^^ = 1 mV. This specific value for the pump 
voltage is sufficient to induce an excitation of the adsorbate spin. If the excitation has to take place 
with the help of intermediate states, e. g. if the transition matrix element does not allow for the 
direct pathway |xi) — ^ IX2), the pump voltage has to be chosen higher according to the energy of 
these intermediate states. We note that in our case the sum of Vpump and K,j.obe ^^ much smaller 

11 






1.0 
0.8 
0.6 

0.4 

0.6 



^^ 0.4 
0.2 



0.2 

^m 0.15 

0.1 

0.05 



CO, 



0.2 
0.15 



1^ 0.1 
0.05 



^, 



0.02 



D-i 0.015 
0.01 
0.005 



O P 




- Wo = 3 eV 

- Wo — 5 eV 

- wq = 7 eV 
uo - 9 eV 



€|v^mBf|^^A 



10 



15 



20 



25 



30 



35 



40 



V^(mV) 



FIG. 3. The population P^ for all spin states of the surface-deposited structure as a function of the applied 
voltage for various coupling strengths uq. The black dashed lines indicate when the voltage is high enough 
as to excite the surface-deposited spin states obeying the energy condition E^ + E^^ — E^ — Ef — eV > 0. 
The bottom row depicts the voltage-dependent wave functions in orbital representation (for uq = 7 eV), 
evaluated for the highest voltage between two black dashed lines. The color coding is as in figure 2. The 
red arrow indicates the expectation value of the spin. Note that when exceeding the first threshold voltage, 
the spin direction flips, as it is primarily connected to the transition \m = +2) -^ \7n = —2). 

than the voltage V3 needed to excited the third or higher excited states. The tip spin can not flip for 
the transition \xi) — > |x2)j hence we can reduce our considerations to a two-level system (TLS). 

The experimentally accessible quantity is the number of tunneling electron A^ during the time 
interval containing the pump and the probe pulses shifted against each other by a time delay At as 
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FIG. 4. The population Pf in the ground state as a function of the voltage V and the coupling strength uq. 
The voltage is bounded such that the only accessible tunneling channel is the first excited spin state of the 
adatom. the vertical white dashed line indicates the energy difference E'f — Ef. We find that the depop- 
ulation process of the ground state is particularly efficient for uq ~ 7 eV. For stronger or weaker coupling 
strengths, the excitation probability decreases. This can be understood by a simple analytical model (see 
appendix). The prediction of this consideration for the optimal value Uq is shown by the horizontal white 
dashed line. 



displayed in figure 5. 

For large | Ai| and At < (5a), there is no overlap of the pump and probe pulse. The number 
of tunneling electrons is therefore just the sum of the electrons transmitted during probe pulse 
A'probc and pump pulse A'pump separately. This number of electrons will be denotes as A'^o = 



N, 



pump 



^probc- For smaller |At| but still negative At, the pump and the probe pulse have a 



finite overlap, i. e. there is a certain time interval where the effective voltage is Vpump + Vprobe- 
If the current-voltage characteristic is nonlinear, the A^ will be different from Nq because of the 
correlation of probe and pump pulse. In our case, the nonlinearity mainly arises from the coupling 
Hamiltonian i/inf 

For At > the pump and the probe pulses have again no overlap, but the pump pulse depletes 
the population of the ground state. The following probe pulse can only access the ground state 
scattering channel, hence the number of transmitted electrons is lower than Nq. With At increas- 
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FIG. 5. A sketch of the pump-probe arrangement of the vohage pulses. In the case (a) the probe pulse 
advances the pump pulse characterized by At < and the cross correlation of both pulses is zero. On the 
other hand, by varying At > as in (b) one can detemiine the relaxation time for the transition |x2) -^ \xi)- 
The probe and the pump pulse both have a duration of 200 ns with the linear rise/fall intervals of 50 ns. 

ing, there is more time for the surface spin to relax back to the ground state, so A^ approaches A'^o 
again for a large At. 

For a quantitative analysis, we note that all quantities are only parametrically time-depended, 
since the time-dependence of ypump(^) and Vprobc(^) is very slow on the typical time scales of the 
system. This means we can use the results of the steady-state considerations from the section III 
and IV. From the tunneling probability for the different scattering channels T^^-iV) we compute 
the corresponding currents as ji 2 = iehki2/'m)Ti^i2)t(^) with ki = ^y27neV/h'^ and k2 = 



^/2m{eV + Ef^^Ef)/h? for eV > E^—E^ and k2 = otherwise. The total current is expressed 
as the sum j = ji +J2 proportional to the total number of tunneling electrons dN in the given time 
interval dt. As the calculations are based on the quasi-static regime, it is still possible to define 
the conductance G{V) oc dj /dV, which then is parametrically time-dependent, as well. As we 
have pointed out in section IV, the excitation of the state IX2) is favored as long as the voltage 
is sufficiently high. For this reason, the average slope of J2{V) is larger than the average slope 
of ji(y), i. e. the conductance jumps to a higher value if eV exceeds i^f — E^ . The nonlinear 
conductance profile is depicted in figure 6. 

We note the population of the ground state and the first excited state evolve only adiabatically 
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FIG. 6. The conductance of the magnetic tunneling junction for various coupling strengths uq. The value 
of G has been normalized to the situation where uq = 0, in which case the conductance is almost constant. 

with resultant voltage Ipump(^) + H>robe(0, because the excitation |xi) — )■ IX2) occurs instantly on 
the time scale of the voltage pulses. With the surface spin driven into the excited state, it undergoes 
i) a precession (the characteristic time is about 2 ps) and ii) relaxes. Since the computation of N 
involves an integration over the time interval of several hundred nanoseconds, the oscillatory part 
due to the precession does not play a role. We can, therefore, describe the dynamics of the surface 
spin in the case of a falling voltage by using the Bloch equations including damping, but at a zero 
frequency. The time scale for this relaxation process, chosen according to the experiment [8], is 
again very slow for our system. 

For a pump-probe delay At < such that the probe pulse still advances the pump pulse, we 
therefore obtain 

-^ 

with j{V) = eh{kiTi^{V) + k2T2^{V))/m. For At > on the other hand, there is no overlap 
of the two pulses and hence the number of electrons transmitted due to the pump pulse is given 
by A^pump = Jl I '^^ J(Vpump(^))- The probe pulse can only access the ground state |xi)- The 
tunneling probability and equivalently the current are in this case proportional to the population 
Pf, which we obtain by calculating PfiVpuLp + VpJbc) < ^ ^^'^ '^hen computing the increase of 
Pf{t) back to 1 due to the relaxation. The number of tunneling electrons due to the probe pulse 



dtj{Vpump{t) + Vprohc{t)) 



(13) 
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FIG. 7. The relative number of the transmitted electrons with respect to A'^o as a function of the delay Ai 
of pump and probe pulses for various interaction strengths uq. The curve has a peak at At = —200 ns 
which refers to a complete overlap of the probe and the pump pulses. Since the population in the ground 
state is lowered by the excitation of the surface spin, N — Nq becomes negative for At > and approaches 
zero again when the excited state relaxes back into the ground state. Our calculations, thus, confirm that 
threshold phenomena may also arise under different excitations: \x2) is excited directly. For comparison, 
the inset shows the experimental result from Ref. 8 where the same final state is reached by a sequence of 
excitations with Am = 1. 



A^probe = Jl I '^^ Pf {t)ji{Vprohc{'t)) Is, thus, lowcred in comparison to the probe pulse advancing 
the pump pulse. 

Figure 7 shows the total dependence of N on the time delay At. Again the curves in figure 7 
evidence a strong dependence on the coupling of the tunneling electron to the local moment uq. 
The peak at At = —200 ns (corresponding to a complete overlap of the pump and the probe pulses) 
is due to the nonlinearity of current- voltage characteristic, as j(Vpump + ^pj^.^) > j{VpuLp) + 
.?( %robc)- ^^^ magnitude of this effect can be explained by means of the the conductance profile 
(figure 6). For small values of uq, the relative jump of the conductance is small, so j{V) is nearly 
linear, i. e. j{Vpnmp{t) + ^probc(t)) > j(^pump) + j(^probo)- Hence, N only slightly exceeds 
A'^o when both pulses coincide. For this reason, the first peak in figure 7 is small for smaller 
Uq. The maximum becomes more pronounced when approaching the optimal value of uq for the 
depopulation of the ground state. From the conductance profile (figure 6) we conclude that the 



relative jump of G when exceeding the threshold voltage attains the highest value for uq 



Ua. 
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The nonlinearity then plays a dominant role, resulting in an especially strong increase of A^ with 
respect to A^o- For larger coupling strength, the current- voltage dependence resembles again the 
linear case. Note that the averaged conductance is much smaller for uq > Uq than in the case 
Uq < Uq. The minimum at At = is caused by the reduction of the ground state population 
Pf due to the advancing pump pulse. We have already discussed the efficiency of this process 
(figure 4). For values uq ^ u^ or uq 3> Uq, Pf < 1 holds. Hence, A^probc and for this matter A^ 
decrease only by a small amount. The depth of the minimum is maximal for uq = Uq. 

We remark that the relative jump in the conductance profile also decreases for /3 < 1. The 
tunneling channel connected to the excitation |xi) — > IX2) is much less effective for | |)-electrons 
emitted from the tip. Hence, the major contribution for this electrons arises from the adsorbate 
spin remaining in the ground state. Hence, the conductance for voltages eV > -Ef — E^ differs 
only slightly from the value left to this threshold. This effect again shows a linear dependence on 
the tip polarization. The total value of the conductance on the other hand increases with smaller 
/3, which is due to the reduction of the effective potential for 1 1) -electrons with the adsorbate spin 
in the ground state. 

Similarly to the discussion in section IV, the special features with respect to Uq become less 
visible if we reduce /3. Hence, the height of maximum in figure 7 decreases. This can be explained 
in terms of the conduction profile which more and more resembles the linear case. On the other 
hand, the transition of the surface spin to the first excited state effectively blocks the 1 1) -electrons 
from the tips and therefore reduces the current more than in the case (3 = 1. We thus conclude that 
the minimum becomes deeper for a realistic tip polarization. 

VI. CONCLUSIONS 

We analyzed theoretically the dynamics of a surface-deposited localized magnetic moment, 
triggered by a spin-polarized current pulse from an STM tip in the presence of a magnetic field. 
We showed and explained how and why the spin population of the surface-deposited structure can 
be manipulated by applying appropriate voltages and investigated the dependence on the coupling 
strength uq of the tunneling electrons to the localized moments. A simplified analytical model 
allowed us to predict a value of uq for which the excitation probability of the surface spin is 
particularly high. The obtained value is in excellent agreement with the numerical calculations. 

The relaxation was incorporated by using the Bloch equations, as our model does not provide 
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a relaxation mechanism. For future work, this point deserves a special consideration. In this 
respect, the main advantage of the Cu2N-overcoat is the magnetic decoupling of the adatom and 
the substrate, the excitation of magnons is strongly suppressed leading to prolonged relaxation 
times. In general, we suggest the spin-orbit interaction, that transfers the energy to the substrate, 
as the major dissipation channel. 

The above treatment is based on an adiabatic approximation in the sense, that the voltage pulses 
durations are very long on the time scale of the electron and spin dynamics, which is a very well 
justified assumption for the experiments under consideration (figure 1). In future experiments one 
may envisage however to go beyond the relaxation dynamics with the aim to image the precession 
of the adsorbate spin as well. In this case picosecond pulses are required. Such pulses have already 
been utilized for tracing the ultrafast spin dynamics in a proof of principle experiments [25, 26]. 
Altematively, one may consider laser-pulse induced tunneling currents, as in the recent experi- 
ments [27, 28]. Generally, the characteristic times for the tunneling process (some femtoseconds) 
and the spin dynamics (picoseconds) are still very different. Hence, the dynamics driven by the 
voltage pulses with picosecond durations can be adiabatically decomposed into (i) the time evo- 
lution of the electron states, (ii) the electron pulse driving the surface spin, and (iii) an adiabatic 
dependence of the current on the surface-deposited spin dynamics (which images of the spin dy- 
namics). This time scale separation is the key for future experiments and calculations on road to 
full insight into the spatiotemporal evolution of the quantum spin dynamics. 



VII. APPENDIX: ANALYTICAL MODEL 

In section IV that for the parameters of interest the excitation probability exhibits a maximum 
for uo ~ 7 eV. To explain the specific value of the maximum and to gain further insight in the 
underlying physics, we introduce here some reasonable simplifications that allow an analytically 
solvable model. Since the scattering channels are limited to the first two states of the surface- 
deposited magnetic moment and the tip electron spin can not flip within this transition, we can 
reduce our considerations to a two-channel model. The operator f ■ S can be replaced by Sz as 
the electron spin is fixed in the z-direction. For small voltages, the gradient of the potential barrier 
can be neglected. We can, thus, assume a rectangular shape of the barrier. The function w{x) 
describing the spin-interaction area is localized around the Fe atom. Hence we approximate the 
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coupling Hamiltonian by 

i^int = uoC 6 (^^) , (14) 

where xq is a typical length associated with w{x), L is the width of the vacuum barrier and C is 
the remaining coupling matrix that contains the transition amplitudes, i. e. C^p = {Xa\S z\xi3) ■ 

The wave functions V^^(x) are then linear combinations of plane waves ^.^^K^-^) for x < or 
X > L and e^'^^^~^^ for a; G (0,L). Therespective wave vector follows as A;j, = ^y2m{eV + Ef — E^)/h?, 



whereas q^ = \/lm{yB — eV^ — -Ef + E^)/fi?. Since the height of the vacuum barrier Vb is much 



greater than the electron energy, we can approximately set q^ = q = ^/2mVB/h'^- For < x < L 
we write ^j'(x) = a^e'^^ + [5^e~''^ and for x > L we define ^^(x) = a^e^^"^ + (5^e~^^"^. The 
coefficients are linked by the corresponding transfer matrices. Imposing the standard matching 
conditions for the wave-function and its derivative we obtain for the reflection coefficient of the 
first excited spin state in the interval (0, L): 

< >mxoUoCi2 , ^> f^h , K-^q-2mxQUQC 22 



where > and < labels the coefficients to the left and to the right of x = L. We can furthermore 
neglect the term k2/q because of the small voltages. The remaining terms are then real, so that the 
absolute value of /3^ is minimal if the right side in equation (15) is minimal. Hence, af becomes 
small near the optimal value u^. We conclude that /3^ is small for q = 2mxoUQC22/h'^, yielding 
for the optimal value Uq = h^q/ {2niXQC 22)- 
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